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Abstract 

Certain unresolved ambiguities surround pressure determinations for 
incompressible flows, both Navier-Stokes and magnetohydrodynamic. For 
uniform-density fluids with standard Newtonian viscous terms, taking the 
divergence of the equation of motion leaves a Poisson equation for the 
pressure to be solved. But Poisson equations require boundary condi- 
tions. For the case of rectangular periodic boundary conditions, pressures 
determined in this way are unambiguous. But in the presence of "no-slip" 
rigid walls, the equation of motion can be used to infer both Dirichlet and 
Neumann boundary conditions on the pressure P, and thus amounts to 
an over-determination. This has occasionally been recognized as a prob- 
lem, and numerical treatments of wall-bounded shear flows usually have 
built in some relatively ad hoc dynamical recipe for dealing with it, often 
one which appears to "work" satisfactorily. Here we consider a class of 
solenoidal velocity fields which vanish at no-slip walls, have all spatial 
derivatives, but are simple enough that explicit analytical solutions for 
P can be given. Satisfying the two boundary conditions separately gives 
two pressures, a "Neumann pressure" and a "Dirichlet pressure" which 
differ non-trivially at the initial instant, even before any dynamics are 
implemented. We compare the two pressures, and find that in particu- 
lar, they lead to different volume forces near the walls. This suggests a 
reconsideration of no-slip boundary conditions, in which the vanishing of 
the tangential velocity at a no-slip wall is replaced by a local wall-friction 
term in the equation of motion. 

1 Introduction 

It has long been the case that pressure determinations for incompressible flows, 
both Navier-Stokes and magnetohydrodynamic (MHD), are known to be highly 
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non-local. Taking the divergence of the equation of motion 

^+vVv = j^-VP + i.VV (1) 
ot pc 

and using V • v = leaves us with a Poisson equation for the pressure P, which 
is said to function as an equation of state: 

V^P = -V • (v ■ Vv - ) (2) 

pc 

Here, v = v(a;, t) is the fluid velocity field as a function of position and time, B is 
the magnetic field, j = cVxB/47r is the electric current density, c is the speed of 
light, ly is the kinematic viscosity, assumed spatially uniform and constant, and 
P is the pressure normalized to p the mass density, also spatially uniform. 
and (H) are written for MHD. Their Navier-Stokes equivalents can be obtained 
simply by dropping the terms containing B and j . 

If we are to solve d) for P, boundary conditions are required. In the im- 
mediate neighborhood of a stationary "no-slip" boundary, both the terms on 
the left of (|l|) vanish and we are left with the following equation for VP as a 
boundary condition: 

VP = ^.VV + (3) 

pc 

We now focus on the Navier-Stokes case, where the magnetic terms dis- 
appear from (|3|), for simplicity. All the complications of MHD are illustrated 
by this simpler case. It is apparent that must apply to all components of 
VP, and that while the normal component of VP is enough to determine P 
through Neumann boundary conditions, the tangential components of (j^) at 
the wall equally well determine P through Dirichlet boundary conditions. This 
is a problem which some inventive procedures have been proposed to resolve, 
usually by some degree of "pre-processing" or various dynamical recipes which 
seem to lead to approximately no-slip velocity fields after a few time steps (e.g. 



[Gresho 1991 , |Roache 1982| and [Canuto et al. 1988| ). It is not our purpose 



to review or critique these recipes, but rather to focus on a set of velocity fields, 
related to Chandrasekhar-Reid functions jChandrasekhar 1961 , for which 



is explicitly soluble at a level where the Neumann or Dirichlet conditions can 
be exactly implemented. In § 2 , we explore the difference between the two pres- 
sures so arrived at. Then in § 3, we propose a replacement for the long-standing 
practice of demanding that all components of a solenoidal v vanish at material 
walls, in favor of a replacement by a wall friction term for which the above 
mathematical difficulty is no longer present. Of course, similar statements and 
options will apply to all comparable incompressible MHD problems. 

2 Pressure determinations 

Restricting attention at present to the Navier-Stokes case, we consider two- 
dimen- sional, solenoidal, velocity fields obtained from the following stream 
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function: 

jp{x, y) = Ckx cos {kx) [cos {\y) + Ak\ cosh {ky)] (4) 

The hyperbohc cosine term in contributes a potential flow velocity compo- 
nent to V which makes it possible to demand that v obey two boundary con- 



ditions: the vanishing of both components at rigid walls | Chandrasekhar 1961 1. 
The function in (^) is even in x and ?/, but can obviously be converted into an 
odd or mixed one by the appropriate trigonometric substitutions. 

The velocity field v = VV^xe^ has only x and y components and is periodic 
in x, with an arbitrary wavenumber k. Ck\ is a normalizing constant, and the 
constants A and A^x can systematically be found numerically to any desired 
accuracy so that both components of v vanish at symmetrically placed no-slip 
walls at y = a and y = —a. In fact, for given fc, an infinite sequence of such 
pairs of A and Akx can be determined straightforwardly. Thus any such v, 
or superposition thereof, is not only solenoidal, but has both components zero 
at y = ±a, and all spatial derivatives exist. Moreover, the "source" term, or 
V • (v • Vv), from the right hand side of (^), is of a relatively simple nature for 
such a V, since every term in it can be written as a product of exponentials of 
kx, Xy and ky. It is straightforward to find an inhomogeneous solution for P, 
which then is the same for all boundary conditions for a given v of the form 
stated. To this inhomogeneous part of P must be added a solution of Laplace's 
equation. This can be chosen so that the total P may satisfy either the normal 
component of (||) at the walls, or the tangential component of it, but not both. 
The determination involves only simple but tedious algebra. 

We illustrate, in figure]^, an arrow plot of the velocity field given by choosing 
k — 7r/2, A = 2.6424 and Akx = .3499, in units of a = 1. The two pressures re- 
sulting from the satisfaction of the normal and tangential components of can 
best be compared by comparing their respective values of VP, since P itself is 
indeterminate up to an additive constant in both cases. In figure ||, we display, as 
an arrow plot, the difference between the pressure gradients associated with the 
velocity field shown in figure [l]. We have rewritten (|l])-(H) in dimensionless units 
for this purpose, with the kinematic viscosity being replaced by the reciprocal of 
a Reynolds number, which may be defined as Re = ((v^)/ (fc^ -I- A^))^/^/i^. Here, 
the angle brackets refer to the mean of taken over the 2-D box, containing 
one period in the x direction and from y — —a to y = a. The value of Re used to 
construct figure || is i?e = 2293, with the dimensionless version of Ckx = 5000 
in (^). The two pressures are similar but not identical. 

In figure |[ a fractional measure of the difference between the "Neumann 
pressure" Pn and the "Dirichlet pressure" Pd is exhibited as a contour plot of 
the scalar ratio 



(VPd - VP, 



N) 



\2 



There is no absolute significance to the numerical value of this ratio. It initially 
increases with Re approaching a maximum of about 2% near the wall for Re > 
10. It is considered interesting however, that the fractional difference is nearly 
x-independent where it is largest. That occurs formally because the algebra 
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Figure 1: Velocity field: v = V^Axe^ using ip from 
and Akx = .3499 



with k = 7r/2, A = 2.6424 



reveals it to be dominated by a term which varies as cosh (4fc?/) — cos (Akx) in a 
region where ky > 1. 

It is amusing but perhaps not significant to superpose the velocity field from 
(||) with a parabolic plane Poiseuille flow of a larger amplitude. The resulting 
flow field is shown in figure ^, and it bears a striking but perhaps not significant 
similarity to the flow patterns seen in two-dimensional plane Poiseuille flow 



[Jones & Montgomery 1994] when linear stability thresholds are approached. 
The pressure gradient difference for this case will be fractionally smaller than 
in figure since pure parabolic plane Poiseuille fiow is a rare case where the 
two pressures happen to agree, and it quantitatively dominates the pressures 
determined from equation (0) in this example. 



3 Discussion and a possible modification 

An alternative to the no-slip condition is the "Navier" boundary condition 
[ Lamb 1932 |: the slip velocity at the wall surface is taken to be proportional to 
the rate of shear at the wall. This may be expressed AV = Lg'y where AV^ is the 
slip velocity of the fiuid at the wall, 7 is the rate of shear at the wall and Lg is a 
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Figure 2: The difference between pressure gradients: VPn — VPd with Re = 
2293 



constant with the dimen sions length. Molecular dyna mic simulations of Newto- 
nian liquids under shear | Thompson fc Robbins 199C ] have shown this to be the 
case under some circumstances. In fact recent work [ [Thompson fc Troian 1997 | 
has shown that, in cases where the shear rate is large, there is a nonlinear 
relationship between Lg and 7. 

We note that the velocity field shown in figure |l| does not lead to one which 
obeys the Navier boundary condition, after an initial time step, where the fluid 
has been allowed to slip at the wall. If the velocity field determined by (^) is 
advanced in time using (Q) with the "Neumann pressure" , the proportionality 
between the slip velocity and the rate of shear at the wall, after the initial time 
step, varies sinusoidally with x. 

It is to be stressed that we are concerned here only with initial conditions, 
not with circumstances under which initial slip velocities might be coaxed dy- 
namically into vanishing after some time. 

It is difficult to see in what sense the velocity field obtained from (^) might 
be an unacceptable one from the point of view of the Navier-Stokes or MHD 
descriptions. It seems to have all the properties that are thought to be rele- 
vant. The family of functions of the same x-periodicity in can be shown 
to be orthogonal, and is a candidate for a complete set, in which any v might 
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Figure 3: Normalized mean square pressure gradient difference: (S/Pd — 
VPAr)2/((VPAr)2) with Re = 2293. Note that the fractional difference between 
the two values of VP is significant only near the wall. 



be expanded, when supplemented by flux-bearing functions of y alone. The 
mathematical question of which if any velocity fields, which are both solenoidal 
and vanish at the wall, would lead to Neumann and Dirichlet pressures that 
were in agreement with each other, must remain open. Indeed, the question of 
whether there are any, without some degree of "pre-processing," must remain 
open. This is an unsatisfactory situation for fluid mechanics and MHD, in our 
opinion, even if it is a not unfamiliar one. The search for alternatives seems 
mandatory. 

One alternative that may be explored is one that seemed some time ago, in 
a rather different context [Shan & Montgomery 1994a, b], to have worked well 
enough for MHD. Namely, we may think of replacing the requirement of the 
vanishing of the tangential velocity at a rigid wall with a wall friction term, 
added to the right hand side of , of the form 

(6) 



r(x) 



where the coefficient 1/t(x) vanishes in the interior of the fiuid and rises sharply 
to a large positive value near the wall. The region over which it is allowed 
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Figure 4: Velocity field shown in figure |l| plus parabolic plane Poiseuillc flow: 
V = VV-xe^ + 2(y2 _ l)e^ with Ck\ = 1 



to rise should be smaller than the characteristic thickness of any boundary 
layer that it might be intended to resolve, but seems otherwise not particularly 
restrictive. Such a term provides a mechanism for momentum loss to the wall 
and constrains the tangential velocity to small values, but does not force it to 
zero. The Dirichlet boundary condition disappears in favor of a relation that 
permits the time evolution of the tangential components of v, while demanding 
that P be determined solely by the Neumann condition (the normal component 
of (^) only). In a previous MHD application [Shan & Montgomery 1994a, b] 
dealing with rotating MHD fluids, the scheme seemed to perform acceptably 
well, but was not intensively tested or benchmarked sharply against any of the 
better understood Navier-Stokes flows. This comparison seems worthy of future 
attention. 

The work of one of us (D.C.M.) was supported by hospitality in the Fluid 
Dynamics Laboratory at the Eindhoven University of Technology in the Nether- 
lands. A preliminary account o f this work was presented or ally at a meeting of 
the American Physical Society [ Kress fc Montgomery 1999 |. 
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